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Q ■ Abstract 



We study LHC phenomenology of mirage mediation scenario in which anomaly and 
modulus contributions to soft SUSY breaking terms are comparable to each other. 



5_i ■ A Monte Carlo study of mirage mediation, with model parameters a = 1, M( 



o — 

500 GeV, um = 1/2, n# = 1 and tan/3 = 10, is presented. It is shown that 
masses of supersymmetric particles can be measured in a model independent way, 
providing information on SUSY breaking sector. In particular, the mass ratio of 
gluino to the lightest neutralino for the benchmark scenario is determined to be 
1.9 < mg/m^,o < 3.1, well reproducing theoretical input value of m g /m^o ~ 2.5 
which is quite distinctive from the predictions m g / m^o > 6 of other SUSY scenarios 
in which gaugino masses are unified at the GUT scale. The model parameters of 
mirage mediation can be also determined from various kinematic distributions. 



1 Introduction 



Weak scale supersymmetry (SUSY) is one of the most promising candidates of new physics 
beyond the standard model (SM) pp. It provides a solution for gauge hierarchy problem 
and complies with gauge coupling unification. Another nice feature of the SUSY theory 
with R-parity conservation is that the lightest supersymmetric particle (LSP) is a natural 
candidate for the non-baryonic dark matter (DM) in the universe. 

SUSY phenomenology crucially depends on the masses of SUSY particles and their 
properties, which might result from spontaneous SUSY breaking in a hidden sector. The 
SUSY breaking is communicated to visible sector through some messenger interactions 
depending on models. Soft SUSY breaking terms of visible matter fields are then de- 
termined such that different SUSY breaking and mediating mechanism leads to different 
pattern of masses and properties of SUSY particles. 

Recently Kachru et al. (KKLT) has provided a concrete set-up of string compactifi- 
cation, in which all moduli are fixed and Minkowski (or de Sitter) vacuum is achieved 
[2]. The KKLT-type moduli stabilization scenario leads to an interesting pattern of soft 
SUSY breaking terms, to which modulus and anomaly contributions are comparable to 
each other if gravitino mass 777.3/2 ~ 10 TeV [3]. A noticeable feature of the mixed 
modulus-anomaly mediation (a.k.a mirage mediation) is that soft masses are unified at a 
mirage messenger scale [I] 



with a representing the ratio of the anomaly to modulus mediation. The mirage messenger 
scale M m i r is hierarchically lower than Mqut for a positive a ~ 0(1). Such a low (mirage) 
unification scale for soft masses leads to a SUSY mass spectrum which is quite distinctive 
from those in other SUSY breaking scenarios such as mSUGRA, gauge mediation and 
anomaly mediation. Some phenomenological aspects of the mirage mediation have been 
investigated by several authors 0, E] • 

In this paper, we investigate the LHC signatures of the mirage mediation performing 
a Monte Carlo study for a benchmark point in the scenario and show that SUSY particle 
masses can be determined in a model independent way, providing some valuable infor- 
mation on SUSY breaking sector. If the weak scale SUSY is realized in nature, SUSY 
particles would be produced copiously at the LHC [7J E], which is scheduled to start in 
2007. Gluinos and squarks, which are directly produced from the proton-proton collision, 
will decay to the LSP and SM particles in the end, with non-colored SUSY particles as 
intermediate states in general. The precise measurement of the masses of SUSY particles 
might be possible with reconstruction of the cascade decay chains [9], [10] . Kinematic edges 
and thresholds of various invariant mass distributions can be measured experimentally and 
then the SUSY particle masses would be determined in a model independent way. In turn, 
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SUSY breaking mechanism might be reconstructed from the measured SUSY spectrum 
and signatures. 

This paper is organized as follows. The basic features of the mirage mediation are 
briefly described in section 2. The Monte Carlo study for a benchmark point is presented 
in section 3. Section 4 is devoted to the conclusions. 



2 Mirage Mediation 

In KKLT-type moduli stabilization scenario, the light modulus T which determines the SM 
gauge couplings is stabilized by non-perturbative effects and the SUSY-breaking source is 
sequestered from the visible sector. The non-perturbative stabilization of T by modulus 
superpotential results in a suppression of the modulus F component: 



F*_ ^ 

~T ~ HM Pl /m 3/2 ] 



m 3/2 



(2) 



which is comparable to anomaly mediated soft mass of O (7713/2 /47r 2 ) for 7773/2 near the 
TeV scale. The soft terms of visible fields are then determined by the modulus mediation 
and the anomaly mediation if the SUSY breaking brane is sequestered from the visible 
sector. The soft terms of canonically normalized visible fields are given by 
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where A a are gauginos, 0j are the scalar component of visible matter superfields $j and 
Uijk are the canonically normalized Yukawa couplings. For F T /T ~ 7773/2 /47r 2 , the soft 
parameters at energy scale just below M GUT are determined by the modulus-mediated 
and anomaly-mediated contributions which are comparable to each other. One then finds 
the boundary values of gaugino masses, trilinear couplings and sfermion masses at Mqut 
are given by 
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where a represents the anomaly to modulus mediation ratio: 

"73/2 
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M ln(M P1 /777 3 /2) 
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with M the pure modulus mediated gaugino mass, while a, and q parameterize the 
pattern of the pure modulus mediated soft masses. The one-loop beta function coefficient 
b a , the anomalous dimension 7^ and its derivative % and 9i are defined by 
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where the quadratic Casimir C%(fa) = (N 2 — 1)/2N for a fundamental representation fa 
of the gauge group SU (N), C%(fa) = q 2 for the U(l) charge of fa, and Wy = Y,m VikiVjki 
is assumed to be diagonal. The explicit expressions of b a , 7$, ji and 9 in the MSSM are 
given in appendix A. In this prescription, generic mirage mediation is parameterized by 



a, M , dj, q, tan/3, 



(7) 



where tan/3 is the ratio of the vacuum expectation values of the two neutral Higgs fields. 

An interesting feature called mirage mediation arises from the soft masses of Eq. (JJ|) 
at Mquti due to the correlation between the anomaly mediation and the RG evolution of 
soft parameters. The low energy gaugino masses are given by jl] 
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implying that the gaugino masses are unified at M m j r , while the gauge couplings are 
unified at Mqut- If the is small or a* + cij + afc = q + Cj + Ck = 1, the low energy 
values of A^-fe and m 2 are given by [I] 
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where Yi is the t/(l)y charge of fa. Therefore, the first- and second-generation sfermion 
masses are also unified at the M mir scale if the modulus-mediated squark and slepton 
masses have a common value, i.e. Cq = c~ e (= cm). 
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Phenomenology of mirage mediation is quite sensitive to the anomaly to modulus 
mediation ratio a as well as the parameters dj and q [H, El [6] . When a increases from zero 
to a positive value of order unity, the nature of the neutralino LSP is changed from bino- 
like to Higgsino-like via a bino-Higgsino mixing region. This feature can be understood 
from the dependence of the gaugino masses on a. The gluino mass M 3 decreases as a 
increases, while the bino mass Mi increases. Smaller M3 then leads to smaller | and 
Higgsino mass parameter at the weak scale through smaller stop mass square. While 
the gaugino masses are not sensitive to a« and q, the low energy squark, slepton and 
Higgs masses depend on those parameters through their boundary values at Mqut and 
their RG evolutions and the mass mixing induced by the low energy A parameters. 

For the original KKLT compactification of type IIB string theory [2], one finds that 
a — 1, Oj = Cj = 1 — Hi, where rij are modular weights of the visible sector matter fields 
depending on the origin of the matter fields. The corresponding mirage messenger scale 
for a = 1 is given by M m j r ~ 3 x 10 9 GeV. Such an intermediate (mirage) unification scale 
which is hierarchically smaller than Mqut, leads to quite a degenerated sparticle mass 
spectrum at EW scale, compared to that for mSUGRA-type pure modulus mediation 
(a = 0). In this paper, we will consider the intermediate scale mirage mediation (a = 1) 
as a benchmark scenario for a detailed Monte Carlo study of LHC phenomenology. 

3 Collider signatures 
3.1 A benchmark point 

We perform a Monte Carlo study for LHC signatures of mirage mediation with the fol- 
lowing model parameters; 

a = 1, M = 500 GeV, a M = c M = 1/2, a H = c H = 0, tan/3 = 10, (10) 

where cm is a common parameter which parameterize the pattern of the pure modulus 
mediated masses for squarks and sleptons and ch for soft Higgses. This choice of model 
parameters is denoted as a blue dot on the (a, Mo) plane of Fig.l, which was taken from 
the figure 12 (c) in ref. [6]. In the Fig.l, the magenta stripe corresponds to the parameter 
region giving a thermal relic density consistent with the recent WMAP observation pTj . 
i.e. 0.085 < QoAih 2 < 0.119. The benchmark point is also consistent with constraints on 
particle spectra and b — > 57 branching ratio. 

The SUSY particle mass spectrum at the electroweak(EW) scale was computed by 
solving the RG equations with the model parameter set (TlOl) . For the gluino and the first 
two generation squark masses, we find 

m § = 884.4 GeV, = 776.0 (771.9) GeV, = 733.5 (741.8) GeV, (11) 

while the masses of the third generation squarks are 

m ~ h 2 = 703.9 (734.6) GeV, = 545.3 (782.0) GeV (12) 
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Figure 1: Parameter space (a, Mo) with an = ch = 0, om = c-m — 1/2 and tan/5 = 10. 
The blue point corresponds to our benchmark point, which satisfy relic density bound 
of WMAP observation and is consistent with other experimental constraints on particle 
masses and b — > 57 branching ratio. 

On the other hand, the slepton and sneutrino masses are 

m &R(L) = 382.0 (431.5) GeV, m fl(2) = 378.9 (435.4) GeV, = 424.1 GeV. (13) 

For the neutralino and chargino masses, we have 

m _o 2 3 4 = {355.1, 416.1, 478.7, 535.6} GeV, m-± 2 = {408.2, 533.5} GeV, (14) 

which correspond to the following bino,wino and Higgsino masses at electroweak scale, 

Mi = 367 GeV, M 2 = 461 GeV, // = 475 GeV (for tan/5 = 10). (15) 

Finally, Higgs masses are given by 

m h = 115 GeV, m H = 528.6 GeV, m A = 528.3 GeV, m H ± = 534.4 GeV. (16) 

Fig. |2](a) shows the mass spectrum for the benchmark point. For comparison, we also 
show the spectrum of a mSUGRA point, which gives the same mass of neutralino LSP as 
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Figure 2: The mass spectrum for (a) the benchmark point in mirage mediation and (b) 
a mSUGRA point. Both cases give the same mass for the lightest neutralino. 

the benchmark point, in Figj2](b). One can notice that the spectrum of the benchmark 
point is quite degenerated, compared to the mSUGRA point. 

With these mass parameters for the benchmark point, the neutralino LSP turns out 
to be bino-like, but with non-negligible wino and Higgsino components. Neutralino pair 
annihilation into a gauge boson pair {ZZ or I/F+IF - ) is then efficient so that thermal 
relic density of the neutralino LSP satisfies WMAP bound though the LSP is rather 
heavy {m^o ~ 355 GeV). 

For the benchmark point, the ratio of gaugino masses at EW scale is given by Mi : 
M 2 : M 3 = 367 : 461 : 850 ~ 1 : 1.26 : 2.32, which is quite degenerated comparing to 
the typical ratio M\ : M 2 : M 3 ~ 1 : 2 : 6 for mSUGRA-like scenarios in which gaugino 
masses are unified at the GUT scale. For the benchmark point, the mass ratio of the 
gluino to the lightest neutralino is given by rrig/m^o ~ 2.5, which is different from the 
typical mSUGRA predictions m~ g / m^o > 6 and therefore implies that the gaugino masses 
are NOT unified at the GUT scale. In the next subsection, we will see the mass ratio 
rrig/m^o can be measured experimentally so that we can get some information on SUSY 
breaking sector. 
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3.2 Monte Carlo Events 



A Monte Carlo event sample of the SUSY signals for proton-proton collision at an energy 
of 14 TeV has been generated by PYTHIA 6.4 [12]. The event sample corresponds to 30 
fb^ 1 of integrated luminosity, which is expected with the 3 year running of the LHC at 
low luminosity. We have also generated SM background events i.e., ti events equivalent 
to 30 fb' 1 of integrated luminosity and also W/Z + jet, WW/WZ/ZZ and QCD events, 
with less equivalent luminosity, in five logarithmic pr bins for 50 GeV < pt < 4000 GeV. 
The generated events have been further processed with a modified version of the fast 
detector simulation program PGS (Pretty Good Simulation) [13] , which approximate an 
ATLAS- or CMS-like detector with reasonable efficiencies and fake rates. 

The total production cross section for SUSY events at the LHC is ~6.1 pb for the 
benchmark point, corresponding to ~ 1.8 x 10 5 events with 30 fb' 1 luminosity. Squarks 
and gluinos are expected to be copiously produced at the LHC, if nig, rrig < 1 TeV. With 
the masses of gluino and squarks in the Eq. ffTTj) . the production cross sections for the 
gg, gq, and qq pairs are about 0.3 pb, 2.7 pb and 2.0 pb, respectively, such that the 
gluino-squark pair production and the squark pair production dominate. 

The produced squarks and gluinos decay generally in multistep. In the recent past, 
the following cascade decay chain of squark has been exploited in detail [91 [TO]; 

(17) 



Ql 



especially in mSUGRA framework. It was shown that it would be possible to reconstruct 
both upper edges for the / + /~, l + l~q, and / ± g mass distribution and a lower edge for the 
l + l~q mass coming from backwards decays of the x° m ^ ne Ql rest frame. Those edge 
values of the mass distributions are given by the following analytic formulae, in terms of 
the particle masses involved in the decay chain; 
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We can also measure the upper edge of the distribution of the smaller of the two possible 
Iq masses formed by combining two leptons with a quark jet [T] 



M™ x (low) 
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for 2m? fl - (m| ; + m|o) > 0. (22) 



From the kinematic edge measurements, the SUSY particle masses might be then de- 
termined without relying on a model. For the benchmark point, we can notice that 
m^o > m e - fl . Therefore, the cascade decay chain (|17[) is indeed open so that the well- 
established method of the kinematic edge measurements can be applied to our case. 

In order to reduce the SM background to a negligible level, we apply the following 
event selection cuts; 

(1) At least three jets with P T1 > 200 GeV and P r2 , 3 > 50 GeV. 

(2) Missing transverse energy, E™ tss > 200 GeV. 

(3) E™ ss /M eff > 0.2, where M eff = P T1 + P T2 + P T3 + P T4 + + £ Zeptons P Tl 

(4) At least two isolated leptons of opposite charge, with Pt > 10 GeV and \rj\ < 2.5. 

(5) Transverse sphericity St > 0.1. 

(6) No b-jets. 

Following the analysis in Ref.[5], we have then calculated various invariant mass dis- 
tributions for the benchmark point. The dilepton invariant mass distribution is shown 
in the Fig. [31 Here, the e + e~ + fi + fi~ — e ± fi T combination was used in order to cancel 
contributions for two independent decays and reduce combinatorial background. The SM 
backgrounds which denoted as blue histogram on the plot, are negligible with the above 
event selection cuts, as expected from previous studies [HIS]. We find a clear end point 
in the dilepton mass distribution. A Gaussian-smeared triangular fit to the distribution 
gives M^ ax = 60.61 ± 0.13 GeV, which is ~ 0.3% lower than the calculated value of 

M rnax = gQ g Ge y for thg decays ^0 _^ J r1 ± . 

The dilepton plus jet invariant mass Mu q has been also calculated by combining the 
dilepton with one of the two hardest Pt jets. The two hardest jets are expected to 
come from squark decays q — > \q as dominant production processes result in a pair of 
squark . The distribution for the smaller of the two possible llq masses is shown in 
Fig. H|(a). The upper edge value of the l + l~q masses is given by M™^ x = 341.5 ± 3.6 
GeV, which is obtained from a Gaussian smeared fit plus a linear background. The 
fitted value is consistent with the calculated value of M^ x = 341.4 GeV for the decays 

QL Xlq -> l%l T q -> x\l + l~q (with m~ qL = 776 GeV). 

Now, we further require that one Mu q should be less than 350 GeV and the other 
greater. Two / ± g masses are then calculated using the combination of l + l~q with the 
smaller Mu q (< 350 GeV). Fig. H^b) shows the mass distribution of the larger of the 
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Figure 3: Dilepton (e + e —e ± fi T ) invariant mass distributions for the benchmark 

point. Blue histogram corresponds to SM background. 

two / =t g system. A fit to the mass distribution near the end point is also shown, which 
yield an upper edge value MJ^ ax (high) = 258.6 ± 4.8 GeV, which is consistent with the 
calculated value of M™ ax (high) = 259.7 GeV. 

The mass distribution of the smaller of the two M\ q is shown in Fig. H](c). The upper 
edge value of the distribution from a linear fit is given by M^ ax (low) = 227. 6 ± 3.3 GeV, 
which is consistent with the calculated value of M]^ ax (low) = 226.5 GeV. 

The lower edge for l + l~q mass can be reconstructed from the larger of the two possible 
l + l~q masses, with additional requirement M™ ax /\/2 < Mn < M^ ax . The resulting 
distribution is shown in Fig. H](d). A fit to the distribution gives M™™ = 145.3 ± 2.8 
GeV, which is consistent with the generated value M™™ = 145.4 GeV. 

From the above five edge measurements for the q^ cascade decays (1T7I) . we can deter- 
mine the masses of q^, Ir, x\ an d X? using the eqs. (I18H22I) . 

With given masses of x? and x\i we can § e t information for gluino and squark masses 
from the following gluino decay into right-handed squark plus quark; 

9 -> QrQ -> XiQQ- ( 23 ) 
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Figure 4: Distributions for (a) the smaller of the two Mu q , (b) the larger of two M\ q , (c) 
the smaller of the two Mi q , and (d) the larger of the two Mu q with Mu > MJ^ ax / 'a/2. 

The maximal value of the two jet invariant mass for the gluino decays (1231) is given by 

<?<? 

At the LHC, gluino would be produced mainly from gq or gg pair production. In order 
to obtain event sample for the gluino decay (1231) . the following event selection cuts are 
imposed; 

(1) At least 3 jets with P T1 > 200 GeV, P T2 > 150 GeV and P T3 > 100 GeV. 

(2) EJp iss > 350 GeV and Ef iss /M eff > 0.25. 

(3) Transverse sphericity St > 0.15. 
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Figure 5: Dijet invariant mass distribution (a) without dilepton and (b) with dilepton. 
(4) No b-jets and No leptons. 

In the selection cuts, we required rather hard Pt3 > 100 GeV in order to reduce 
qq production events. After the selection cuts, we calculate dijet invariant mass M qq 
using the three hardest jets. The smallest M qq among the three possible combinations is 
then shown in Fig. [5] (a). The M qq distribution is fitted near end point by a Gaussian 
smeared linear function with a linear background. The resulting edge value is given by 
M™ ax = 431.1 ± 32.1 GeV, which is consistent with the calculated one M™ ax = 432.4 
GeV (with m qR = 733.5 GeV). 

Additional information on m g and m qR might be provided by the following cascade 
decay of gluino, 

g ^q~R.q^ xIqq H j&lho), ( 25 ) 

for which the upper edge value of the two jet invariant mass distribution is given by 



JVC ax (ZZ) 



(mj - m 2 qR )(m 2 qR - m 2 -o] 



1/2 



(26) 



In the gluino cascade decay (J25l) . we consider the right-handed squark decay into x° rather 
than xl, where the \2 further undergoes the dileptonic decay. In the benchmark point, 
branching ratio BR(q~fj — > XzQ) = H% is comparable to BR(qL — > x^l) = 18.5%. This 
is because the gaugino masses are quite degenerated (i.e. M 2 /Mi ~ 1.26) at EW scale 



11 



so that X2 has sizable bino-component in the benchmark point, which is in contrast to 
the typical mSUGRA case where x® * s almost wino-like and therefore BR(q~R — > x^q) is 
negligible. Furthermore, the gluino branching ratio of BR(g — > q^q) ~ 19% is larger 
than BR(g —> q L q) ~ 11.2% for q=u,d so that the decay chain of g — > q R q — > x 2 qq is 
comparable to that of g — > q^q — > x^l- Considering the squark masses that rriq R < rriq L , 
the upper edge value of two jet invariant mass for X2QQ events is essentially determined 
by the Eq. ( 1261) . which involves rriq R rather than rriq L . 

In order to select events which include decay chain ( 1251) . we require 

(1) At least 3 jets with P T1 > 100 GeV and P T2>3 > 50 GeV 

(2) EJp iss > 200 GeV and M eff /E™ iss > 0.2 

(3) At least two isolated leptons with opposite charge 

(4) Mu < 61 GeV and the smallest M Uq < 350 GeV 

(5) Transverse sphericity St > 0.1 

(6) No b-jets 

After the selection cut, dijet invariant mass was calculated using the jet, which gives 
the smallest Mu q , with third and fourth (if any) energetic jets. The smaller of the two 
possible dijet mass is then shown in Fig. El(b). A fit to the distribution gives M™ ax (ll) = 
405.7 ± 8.6 GeV, which is consistent with the calculated value of M™ ax (ll) = 406.9 GeV 
(with rriq R = 733.5 GeV). 

The upper edge measurements of two dijet invariant masses provide lower limits on 
the gluino and squark masses and a strong correlation between rrig and rriq R , but without 
constraint on the upper limit on the masses. In order to further constrain the right- 
handed squark mass, we consider squark pair production and their subsequent decay into 
two quarks plus two LSPs; 

qm&R -> qiXiQ2Xi, (27) 
and construct another variable Mt2 [IS], which is defined by 

M T2 = x mi n m . ss [max{m|(p^\p*J,7?4(p^,p£ 2 )}] , (28) 

where 

m 2 T2 (p q T , k*) =m 2 q + m 2 x + 2(E T E* - p q T ■ k*), (29) 

E q T = \J\Vt\ 2 + ™1 EX = y/\\4\ 2 + m? x . (30) 

Here, p^ 1 and p^ denote the transverse momentum of quark-jets from the squark decays 
and p™ ss is the observed missing transverse momentum and m x is an estimate of the 
lightest neutralino mass. The M T2 distribution has an endpoint at the squark mass rriq R 
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Figure 6: (a) The My 2 distribution with the input value ofm x = 350 GeV as an example, 
(b) A relation between rriq R and m^o from the Mt2 analysis. 

when the input m x is equal to the correct m x o value. In general we obtain a relation 
between and rn x o. 

In order to have event sample for the squark pair production and their subsequent 
decay fl27j) . the following event selection cuts are required; 

(1) At least two jets with P T1 > 300 GeV, P T2 > 50 GeV 

(2) Ef ss > 200 GeV, M eff /E^ iss > 0.3 

(3) Transverse sphericity St > 0.15 

(4) No leptons, no b-jets. 

Fig. [6](a) shows an example of Mti distribution with an input value of m x = 350 
GeV. From a fit to the distribution, we obtain = 749 ±12 GeV for the example. A 
general relation between and m x o from the M?2 analysis is shown in Fig. Mjo). 

Now we can determine six SUSY particle masses, m x o, m x o, rnj R , m^, nig and m,- fl , 
from the various kinematic distributions we considered so far. In order to scan possi- 
ble values of SUSY particle masses, random numbers for m x o, rnj R , nig L , rrig R , and rrig 
were generated within some ranges around their nominal values, while m x o values were 
calculated with the measured MJ^ ax value. The chi-square from the various kinematic 
observables with their errors was calculated to determine the probability for each set of 
masses. Fig. [7] shows the probability distribution of (a) m^o, (b) m x o —m x o, (c) — m^o 
(d) rriq L — m x o, (e) — m x o and (f) m g — m x o, respectively. While the xl mass is de- 
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Figure 7: Probability distribution of (a) m^o, (b) m^o—m^o, (cjm^—m^o, (d)mg L —m^o, 
(e) rriq R — m^o and (f) m g — m^o. 



termined as 

m -o = 3561^° GeV, (31) 
the differences between sparticle masses are rather well constrained, 

m,o - m*o = 60.8 ± 0.1 GeV, mr - m,o = 27.0 ± 0.6 GeV, m Sr - m,o = 436 ± 22 GeV, 

X2 X\ ' Ir Xi ' 1L X\ 

niq R — m^o = 393 ± 16 GeV, m g — m^o = 536 ± 15 GeV, (32) 

respectively. The central values for the estimated masses are consistent with the generated 
ones. And the range for the mass ratio between gluino and the lightest neutralino is given 
by 

1.9 < —5- < 3.1 {reduced \ 2 < 1), (33) 
14 



which is quite distinctive from the typical predictions i.e., m g /m^o > 6 of the other SUSY 
scenarios in which gaugino masses are unified at GUT scale. 
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Figure 8: Probability distribution of (a) a, (b) Mq, (c) cm, and (d) ch- 

So far we have considered 'model-independent' measurements of SUSY particle masses. 
We can also determine the model parameters of mirage mediation scenario from the 
kinematic endpoint measurements. As mentioned in section 2, generic mirage mediation is 
parameterized by M , a, a^, Cj and tan/3. Here, we assume that visible sfermion fields have a 
universal parameter cm(= «m) while the corresponding parameter for Higgs fields is given 
by ch(— clh)- The model is then described by five parameters M ,a,CM,cn and tan/3. 
From the Eq. ([8]) and (Q we can see that the gaugino and the first generation sfermion 
masses at low energy scale are essentially determined by Mo, a and cm- Therefore, we 
can determine the model parameters Mo, a and cm from the measured m g ,mq and mj . 
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With given values of Mo, a and cm, the parameters ch and tan/3 might be obtained from 
the measured m^o and m^o. We have scanned five model parameters in certain ranges 
and calculated the x 2 from the calculated and measured kinematic edge values, in order 
to constraint the model parameters. Fig. [H] shows probability distributions for «,M ,Cm 
and Ch, resulting in the following values of the parameters: 

a = 0.98 ± 0.03, M = 492 ± 22 GeV, c M = 0.48 ± 0.01, c H = 0.02 ± 0.05. (34) 

The model parameters are determined quite accurately, except tan/3 whose value is not 
so well constrained because neutralino masses depend on it rather mildly. The resulting 
central values of the parameters well agree with the input values of Eq (flOl) . 

4 Conclusions 

In this paper we have investigated LHC signatures of mirage mediation by performing a 
Monte Carlo study for a benchmark point in which anomaly and modulus contribution 
to soft SUSY breaking terms are comparable to each other. 

The benchmark point has characteristic features of mirage mediation such as rather 
degenerated gaugino masses at EW scale (Mi : M2 : M3 ~ 1 : 1.26 : 2.3) and en- 
hanced Wino- and Higgsino- components of the lightest neutralino, comparing to typical 
mSUGRA case. Non-negligible wino and Higgsino components of the LSP leads to right 
amount of thermal relic density which is compatible to WMAP data, though the LSP 
mass is rather large (fn x o ~ 355 GeV). Rather degenerated gaugino masses at EW scale 
lead to the mass ratio of gluino to the lightest neutralino rrig/m^o ~ 2.5, which is quite 
distinctive from the typical predictions i.e., m g /m^o > 6 of the SUSY scenarios which 
have gaugino mass unification at the GUT scale. 

For the benchmark point, the 'golden' decay chain of squark (fT7|) is open, allow- 
ing model independent measurement of SUSY particle masses such as m,^, m ; ~ fl , m^o and 
m^o. We also determined gluino and squark masses from various gluino and squark 
decays. The SUSY particle masses are then determined in a model independent way, 
giving valuable information on SUSY breaking sector. In particular, the measured ratio 
of 1.9 < nig/m^o < 3.1 well reproduce the theoretical input value m^/m^o ~ 2.5 of the 
benchmark point. Therefore, the benchmark scenario may be distinguishable experimen- 
tally from other SUSY scenario in which gaugino masses are unified at the GUT scale. 
Model parameters were also obtained with small errors, from a global fit to the observ- 
ables with their estimated errors. The resulting central values of the model parameters 
well agree with the input values for the benchmark point. 
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Appendix A. 

The one-loop beta function coefficients b a and anomalous dimension 7, in the MSSM are 



given by 



-3 



1 



b 



33 



= 5' 




(35) 



where g 2 and g Y = y/S/Egi denote the SU(2)l and U(l)y gauge couplings. 
The 9i and 7$ which determine the soft scalar masses at Mqut are given by 

®h u = 3gi + g Y -6y?(a Hu + a Q3 + a U3 ), 

®H d = 3#2 +g Y -6yl(a Hd + a Qa + a Ds ) - 2y 2 T (a Hd + a La + a Ea ) 

0g a = yfll + 3 ^ + q9y ~ 2 (y 2 ti a H u + a Qa +a U3 ) + y 2 b {a Hd + a Qz + a D3 )j5; 

. 16 16 2 a 1l \r 

Qua = y#s + y 9y ~ Ayt ™ + a ^ + a us) 6 ^ 

®D a = y#3 + ^9y ~ 4yt(a Hd + a Qa + a D3 )S 3a , 

= 3gj + g Y - 2yl(a Hd + a Ls + a E3 )5 3a , 
°E a = 4g Y - ^y 2 T {a Hd + a L[1 + a E3 )S 3a , 



(36) 
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and 



3 4 . 11 4 o 2j 

7^ u = t; 92 + Y 9y ~ Vt Vt ' 
3 „ 11 



3d • 



7H d = ^2 + y9y- 3 %^-2/A,. 
3 11 

88 

7c/ a = -8^3 + —g Y - ZytKAa, 
22 

7D a = -8#3 + y^y - 2 Vb b y b 5 3a, 

3 11 

7L a = + y#y - vlKha, 

j Ea = 22g*-2y 2 T b yT 5 3a , (37) 



where 



7 1^ 2 o 2 2 2 2 

&w = -y#3 - 3 ^2 - y£y + %t + 

K = -y^-3^ - gflr + Vt+ 6 Vb +y 2 r, 

K = -3g 2 2 -3g 2 Y + 3y 2 b +Ay 2 T . (38) 
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